Nearly Localised States in Weakly Disordered Conductors. 
II. Beyond Diffusion Approximation. 
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We use optimal fluctuation method for a new ballistic cr-model to study the long time dispersion of 
conductance G(t) of a mesoscopic sample. In the long time limit the conductance of a d-dimensional 
sample decays as exp (— A\n d t). At shorter times the new results match those in our previous 
paper Jl). It is found that at very long times the diffraction effects are important and the ballistic 
treatment is not valid. We also suggest a physical picture of trapping. 



I. INTRODUCTION 

This article is a continuation of our paper |TJ where 
we suggested an optimal fluctuation method to study ex- 
ponentially rare fluctuations in weakly disordered con- 
ductors. The quantity of interest was expressed using 
the super-matrix nonlinear cr-model and the result- 
ing functional integral was evaluated by the saddle point 
method. In this way we obtained an intermediate asymp- 
tote of conductance time dispersion in one- and two- di- 
mensional metals. It was, however, pointed out that this 
method has only a limited scope of validity, e.g., it fails 
to describe satisfactorily three-dimensional conductors. 
The failure occurs when the super-matrix corresponding 
to the optimal fluctuation changes rapidly over the mean 
free path I. Then the diffusion approximation breaks 
down and the use of the nonlinear cr-model m can no 
longer be justified. 

In this article we treat the problem using the recently 
proposed ballistic a- model jjj: the generalisation of the 
standard one that correctly accounts for large gradients. 
To be specific, we consider the long time asymptote of 
conductance dispersion. If a voltage V(t) is applied to a 
sample, the current through it is given by the Ohm law 



I(t) = / G(t - t')V(t')dt' , 



(1) 



and we are interested in the behaviour of the conductance 
Git) in the long time limit. 

We first discuss this problem under conditions when 
the diffusion approximation is valid. The super-matrix 
theory is defined by the the partition function (see 
and y for review): 

Z = J T>Qe~ F , F=^str J dr{D{VQf + 2iujAQ}, 

(2) 

where the functional integral over super-matrices Q is 
subjected to the constraint 



The expression for the averaged conductance G(t) looks 
as follows 

G(t) = G e-* /T + / ^e"*-* / VQP{Q}er F (4) 



In Eqs. (g), (Q) v is the density of states, D is diffusion 
coefficient and r is the mean free time. The strategy sug- 
gested in Ref jij consists in studying the condition of the 
extremum of the free energy F in Eq. (0) 



2L>V(QVQ) + iu[A,Q] = 



(5) 



together with the condition at the boundary T between 
the mesoscopic sample and a bulk electrode 



Q\r = A, 

and a self-consistency condition 
4tA f dr 



■kTi 



V 



str{AQ}; A 



(6) 



(7) 



that arises after integrating out u> in Eq. (0). After solv- 
ing equation (^) with boundary conditions ^) , expressing 
frequency u> through time t, using self-consistency con- 
dition ([?]) and substituting the solution Q(r) to the free 
energy (H), we obtain the conductance G(t) with expo- 
nential accuracy. 

Long time retardation in electric response G(t) is 
caused by relatively improbable quasi-localised states 
that are weakly coupled to the bulk electrodes and have 
life time of the order of t. Since the mean square of the 
wave function |\l/(r)| 2 is connected to the super- matrix 
Q(r) via 



|*(r) 



str{AQ}, 



(8) 



= 1. 



(3) 



the self-consistency condition (^) can be regarded as a 
relation between the life-time of a quasi-stationary state 
and its wave function. The super-matrix Q is fixed at 
the boundary of the sample (see Eq. ^)), so, to satisfy 
the self-consistency condition (^) str{AQ} must grow to- 
wards the middle of the sample. The gradients of the 
Q-matrix increase simultaneously with the delay time t. 
Since theory (0) correctly accounts only for the lowest 
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term in VQ, it cannot be used for sufficiently long times 
t. 

The importance of high gradients for long time asymp- 
totes and tails of distribution functions was first an- 
nounced by Altshuler, Kravtsov and Lerner (AKL) ||, 
who found the growth of the corresponding invariant 
charges under renormalization group flow. The optimal 
fluctuation method has been used recently by Falko and 
Efetov [|| and by Mirlin who studied the influence 
of the nearly localised states on the distribution of wave 
function amplitudes || and local density of states |7j. 
These authors found that the gradients of Q become large 
near the centre of the sample and, therefore, the cr-model 
description breaks down. 

To treat consistently this problem, along with the oth- 
ers where ballistic motion of electrons is essential,we sug- 
gested a new version of the nonlinear a- model [|3|. This 
theory operates with a super-matrix distribution func- 
tion 5n(T*), where n is the unit vector of the electron 
momentum direction (p = npp)- It effectively accounts 
for the infinite series in IV Q of which only the leading 
term is kept in the action (|J), and, therefore, correctly 
describes fluctuations of the super-matrix Q with wave 
vectors q ~ l/l. The theory is still restricted, however, 
to the region of validity of semi-classical approximation. 
The condition of extremum for the new action is related 
to the kinetic equation in the same manner as Eq. (j^) 
is related to the diffusion equation. We solve the kinetic 
equation for one- two- and three-dimensional geometries 
and obtain very long time asymptotes of conductance. 

The results are summarised in table Q. The exponen- 
tial decay law for times to <C t <C h/ A can be obtained 
analysing the time dependence of weak localisation cor- 
rections. The "ballistic" part of the long time asymp- 
tote at t > tf, was first studied by AKL who used the 
frequency representation and expanded the conductance 
G(lo) in powers of w. They estimated the growth rate 
of coefficients in this expansion and pointed out the im- 
portance of high gradients. Unfortunately, the Fourier 
transformation to i-representation in two-dimensions was 
carried out with insufficient precision (see for discus- 
sion) and the factor g was lost under the sign of In 2 . 
The AKL approach also failed to predict the existence 
of intermediate asymptote at ?iA _1 <C t <C i& which was 
discovered by the authors |l| using the optimal fluctu- 
ations method for the diffusive cr-model. We found the 
region of validity for this intermediate asymptote and ob- 
tained some estimates for even longer times. These esti- 
mates were recently improved by Mirlin Q who imposed 
somewhat arbitrary effective boundary conditions on the 
super-matrix Q at the point where the diffusion approxi- 
mation breaks down. The value of the action turned out 
to be rather insensitive to the exact form of these bound- 
ary conditions which enabled Mirlin to rederive the AKL 
result in 2d and obtain exp(— A\n 3 1) asymptote in 3d, 
although without the value of the coefficient A in the 
exponent. 

The full ballistic treatment presented in this article 



gives the coefficient in the exponent for 3d; confirms the 
AKL result in 2d; and discovers a new regime in the case 
of a thick wire. We have also found a restriction on the 
validity of ballistic treatment. It turns out, that at ultra- 
long delay times the super-matrix distribution function 
depends strongly on the direction n of the momentum 
developing sharp features with characteristic width Scj). 
At times t ~ £q this width becomes comparable with the 
diffraction angle 5<j) q ~ A/a, where a is a typical size. 
At longer times the diffraction effects become important 
and the ballistic treatment is no longer valid. The values 
of times tQ are presented in Table |[ Apart from deriva- 
tion of these results from the first principles the paper 
describes the optimal fluctuation of random potential in 
a strictly one-dimensional wire that traps an electron for 
time t. We believe that the same mechanism is responsi- 
ble for long time delays in higher dimensions. 

The material is organised as follows. In section II 
we present physical motivations behind the ballistic cr- 
model, find its condition of extremum; and show how 
the ballistic description transforms into the diffusive one 
when the gradients are small. We also discuss the geom- 
etry of the super-matrix distribution function g n (r) and 
introduce its convenient parametrisation. The optimal 
fluctuation method is described in section III. In section 
IV the solution of kinetic saddle-point equation is found 
and the long time asymptote of the conductance is eval- 
uated in two and three dimensions. The one dimensional 
case is discussed in section V. Physical picture of trap- 
ping is presented in section VI. Finally, in section VII we 
discuss the results and the limits of their validity. 



II. EFFECTIVE ACTION FOR QUANTUM 
BALLISTICS 

A. Outline of Derivation 



present a generalised non-linear 
3L which is valid in the ballistic 



In this section we 
super-matrix cr-model 
regime. 

If the quantum effects are neglected, the ballistic 
regime is described by the Boltzmann kinetic equation 
for the distribution function fn(r) of coordinate r and 
momentum p = npp. The quantum description operates 
with density matrix g n (r). To enable averaging over dis- 
order, gnir) should be a super-matrix |3J, analogously to 
the matrix Q in the standard cr-model (||). The quantum 
generalisation of the kinetic equation has the form M 



2vn 



dg n {r) 
dr 



^\^l,fc(r) 



with the additional constraint 



(9) 



(10) 



which is similar to the one imposed on Q. Equation (|^) 
is the required ballistic generalisation of the saddle-point 
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equation (0) and serves as an extremum condition for 
the ballistic action we are constructing. To generate the 
first derivative in Eq. (0) , the action has to have a Wess- 
Zumino type term 



Mg n } = 



du str ( g n (r, u) 



g n {r,u = 0) = A; 



du ' dr 
(r,u= 1) = g n (r), 



dr (11) 
(12) 



where an arbitrary smooth interpolation can be chosen 
as g n (r,u) p0| , and the angular brackets denote aver- 
aging over directions of n. The functional derivative 
SW/5g n (r) is taken with the constraint ( fiol) which guar- 
antees that gnSg n + 5g n g n = and an arbitrary variation 
5g n has the form Sg n = [g ni in]- As a result of variation 



SW = 4 dr str 



dgn 
n——a r 
or 



(13) 



the first derivative appears. There is another way of writ- 
ing the functional W which employs the decomposition 
•h, I'M' ': 



W{g n } = 4 J drstr(Af/- 1 n^) 



(14) 



This representation can be verified by comparing the 
variation of W with Eq. (|l3|). The quantum ballistic 
partition function Z can be presented as an integral over 
distribution functions g n (r) with effective action F: 



Z= Jvg n (r)e- F , 



(15a) 



F = ^J drstr|^A(. 9 (r))-^(.g(r)) 2 | 



TTVVp 

8"~ 



W{g n }. 



(15b) 



The details of the derivation can be found in Ref || . 

Field theory (|l5|) enables us to study any chaotic prob- 
lem, for which the semi-classical approach is valid, irre- 
spective of validity of the diffusion approximation. If 
space gradients are small (IV ' g <C g), the standard treat- 
ment recovers (see H) the Q-matrix theory (|J). To show 
this we expand the matrix g n into a sum over angular 
harmonics keeping only the zeroth and first terms: 



Q(r) 1 



QJ 2 



J(r) • n 



The constraint g 2 = 1 now reads 



Q 2 



1. 



QJ + JQ = 0. 



(16) 



(17) 



Substituting Eq. (|T|) into Eqs. and using conditions 
([l?]), we obtain the partition function in the form 



Q.J 



T 



dr str{iuiAQ -i 



Vf 



Gt 



(VQ)QJ} (18) 



The Gaussian integral over J in Eq. (18) is dominated 
by the vicinity of 



J = l(VQ)Q 
and leads finally to Eq. (0). 



(19) 



B. Symmetries of g-matrices 



Both Q and A are 8x8 matrices which act on the 
8-component super- vectors ^> with the basis H : 



* T = (xi,Xi,<Si,<Si,X2,X2>< s '2,<S , 2), 



(20) 



where \ are the fermionic and S are bosonic variables; in- 
dices 1, 2 correspond to the retarded and advanced Green 
functions with energies Ef ± uj/2. For a super-matrix M 
the super-trace is defined as follows: 

str M = M n + M-22 - M 33 - M Ai + 
M 55 + M m - M 77 - M 88 

In this basis the super-matrix A has the form 

A -\A \ — I 1; i = 1 — 4 
Aij-Aidij, A i -|_ 1 ^ . = 5 g 

For our purposes it is more convenient to use another 
basis where each variable is stands next to its charge 
conjugate partner and fermionic variables are separated 
from bosonic: 



* T = (xi,X2,X2,Xi,Si,S2,S 2 ,S1). 



(21) 



In this new basis, which will be used everywhere in this 
paper, the matrix A has the form: 



A = 



cr 3 <£> r 3 
cr 3 ® t 3 



(22) 



where the blocks in Eq. (|22j) correspond to the fermionic 
and bosonic sectors. The matrices r act inside the 2x2 
blocks, while matrices the cr mix these blocks inside the 
4x4 sectors. 

We first discuss the symmetry of the Q-matrix from the 
standard cr-model (g) focusing on the properties of the 
boson-boson (B) and fermion-fermion (F) sectors only. 
We consider oly the unitary ensemble when an additional 
condition 



[Q b ' f ,1®t 3 ] =0 



(23) 



is fulfilled. Together with constraint (Q), requirement 
(p3|) permits to parametrise the matrices Q B,F with four 
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complex 3-vectors l B ' ,m B,F subject to the constraint 



l 2 = m 2 



Q 



1: 

B.F 



1 +T 3 



i Irr 



1-T 3 



Parametrisation ( p4[ ) is still too general because the ma- 
trix Q has additional symmetries: charge neutrality 



Q = CQ T C T 



Q, 



(25a) 



with the charge conjugation matrix 
ioi ®r\ 



C = 







and pseudo-hermiticity 

Qt = KQ+K = Q, C 



-<J2 ® Tl 



1 

cr 3 ® r 3 



Requirements ([25a|) are satisfied when 



Z* = -m B , Z* = -m F , 



and 
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i B = a(0*> i* = (i F y, a = I o -i o 

1 



(25b) 

(26a) 
(26b) 



Thus, the fermionic sector is parametrised by a real vec- 
tor l F subjected to constraint (l F ) 2 = 1, i.e. a sphere 
S 2 



S 2 = I 1 



l 2 + l 2 2 + l 2 = l, ImZ = 



(27) 



while the bosonic sector is represented by the vector Z 
with two imaginary components if and l B and one real 
if . Due to constraint I =1 the bosonic sector is repre- 
sented by a hyperboloid TL 2 ,: 



Hi 



[l B , -11^-11^ + 11 = 1, Z = /tZ*}. (28) 

The matrix g n of the ballistic er-modcl (|l5| ) obeys the 
constraint gl — 1 and, therefore, can be parametrised 
with Eq. (p4|). The condition of charge neutrality ]Tl| ] 
(j25a|) must be replaced by the generalised version 



5V, 



- C^C T = <?_ 



(29) 



because the charge conjugation changes the sign of the 
W-term in action (]l5|). 

The pseudo-hermitian transformation g n — * Kg+K 
not only changes the sign of the W-term but also replaces 
uj by — Cij* . In the long time asymptote calculation the 
frequen cy u is purely imaginary and the generalisation 
of Eq. (25b) reads: 



Kg+K = g. 



(30) 



Symmetries (|29|) and ( |30| ) impose a new restriction on 
the vectors I and m (compare with Eqs. (|26|)): 



(24) ^ = r = -m* n , = A(^)*, Z^ B = (C) 



(31) 



Since conditions (31) are imposed on the componets 
of two vectors I and m, they are less restrictive than 
Eqs. (|25|). Both g B and g F can be parametrised with a 
complex unit vector l n = + 277^: 



1 +T 3 



<rZ T 



1-T 3 



) crZ 



with the constraint 



i 2 = e - v 2 + m-n) = 1. 



(32) 



(33) 



The geometric meaning of Eq. (|3^) can be described 
as follows. The vector £ = is characterised by 
its absolute value £ and a unit vector v, which corre- 
sponds to a point on a sphere S 2 . Due to the condition 
= $,tj = ^vq, the vector r\ belongs to the plane tan- 
gential to the sphere S 2 at the point v. The condition 
£ 2 — T7 2 = 1 means that at every point of the sphere 
there is an upper part of the two sheet hyperboloid 
with the axis along the radius of the sphere. In other 
words, the vector Z subject to constraint (|33|), belongs to 
the fibre bundle with the base S 2 and the fibre H 2 . The 
Q-matrices in this geometric picture are represented by 
the sub-manifold of this fibre bundle: Q F belongs to the 
base and Q B lies on the fibre over the North Pole of S 2 
(see Eq. (p7|)). The A-matrix corresponds to the bottom 
of TL 2 , in the fibre over the North Pole. 



III. STEEPEST DESCENT PROCEDURE FOR 
QUANTUM BALLISTICS 

The long time asymptote of the conductance G(t) is 
found following the procedure outlined in the introduc- 
tion. The generalised version of Eq. (Q) has the form 



(34) 



where F is the ballistic action ( |15b| ), and the explicit 
form of the functional P is irrelevant within exponential 
accuracy. The functional integral ([34j) over g n is evalu- 
ated using the steepest descent method, which consists 
in: 

1. finding a solution g n (r) of the saddle point equation 



2. expressing the frequency lj through the time t using 
the self-consistency condition 
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At A f dr r A , ^ 



(35) 



which arises as a result of integration over u> in 
Eq. (H); 

3. substituting g n (r) in Eq. (|34|) and obtaining the 
result with exponential accuracy. 

We do not specify boundary conditions for Eq. (||) 
bearing in mind that far from the centre of the sample 
the space gradients become small and g n (r) approaches a 
solution of diffusion equation (|^) . It was shown in Ref. Jl]] 
that for the long times t ^S> h/A the solutions of diffusion 
equation (J^j can be written in the form: 



Q = U 



CT3 <g> T 3 

((T3 cosh 9 + i<72 sinh ( 



where the real angle 9 obeys the equation 



-sinh 9 = 0; % ad = 



(36) 



(37) 



and the constant super-matrices U commute with A. 
They do not enter the action and thus will be omitted. 
As can be seen from Eq. (|3^), only the bosonic block Q B 
has a non-trivial dependence on the coordinates. The 
same holds true in the ballistic region, and from now we 
consider only the bosonic sector of the theory. 

Combining the diffusion asymptote (^) with 



Eqs. 



we find that the vector l n that 



parametrises g n approaches the limit 

li = —lri^79, l 2 = ismhO, l 3 

far from the centre of the sample. 
Eq. (S§) that 



coshf? 



(38) 



One can see from 



Im U — Re L = Im L = 



(39) 



It turns out that the condition (|39j) remains valid even 
in the ballistic regime. Therefore, for the solution of ki- 
netic equation (^) the complex vector Z„(r) is restricted 
to the real three dimensional subspace (p%. The addi- 
tional constraint I 2 = 1 thus acquires the form 



l\-\l 2 \ 2 + ll = l 



(40) 



and defines a one sheet hyperboloid. Convenient coordi- 
nates on the subspace (39) are associated with the cone 
'l — l^l 2 + = and described in appendix [a|. 

Due to the symmetry (pl|), l\ is an odd and ?2,3 are 
even functions of n. Averaging Eq. (H) over directions of 
n and using the parametrisation (p2|), we arrive at the 
continuity condition 

2w F div(n/i) + u{l 2 ) = (41) 

which reduces to the conservation law for the current: 



J=(nh), divJ = 0, 
if the frequency is small and 7 = iur — > 0. 



(42) 



IV. SOLUTION OF KINETIC EQUATION AND 
LONG TIME ASYMPTOTE OF CONDUCTANCE 



In Ref I j| we considered the solution of the diffusion 
equation (37J) either for a disordered wire (ID), or for a 
disc (2D), or for a droplet (3D) of radius R (see Fig |l|) . 
In analysing the kinetic equation (|^) we stick to the same 
geometry. In this section we find the solutions of Eq. (||) 
in space dimensions 2 and 3, postponing the discussion 
of one-dimensional case to the next section. 

Long times correspond to small frequencies to (7 <SC 1). 
This means that the terms containing 7 in the kinetic 
equation can be neglected everywhere, except the cen- 
tral part of the sample. It can be seen from Eq. ( |4^ ) 
that the total flux of the current J is conserved in the 
outer area. Since the current J is directed along the 
radius, flux conservation means that in two- and three- 
dimensional samples the current density decays towards 
the outer boundary together with the space gradients of 
all relevant quantities, and the solution of the kinetic 
equation approaches the diffusion asymptote. The situ- 
ation is, however, different in the one-dimensional case 
when the current and gradients do not decay and the 
diffusion regime is reached only outside the sample. 



A. Qualitative description 



In the space dimensions 2 and 3 the qualitative be- 
haviour of the solutions is the same for both the diffu- 
sive and kinetic equation. We illustrate it with the dif- 
fusion Equation (|37]) which can be regarded as describ- 
ing a chemical reaction where 9 is the concentration of 
a reacting agent. The term DV 2 9 describes the propa- 
gation of the agent in a porous medium and the rate of 
the agent reproduction depends on its concentration as 
7 sinh 9. Since 7 <C 1, generation takes place only in the 
central part of the sample, where the value of 9 is high. 
We refer to this central zone as the "zone of reaction". 
In the outer part of the sample the agent diffuses freely 
(the "run-out zone"), i.e. 9 is a solution of the Laplace 
equation V 2 # = 0. An azimuthally symmetric solution 
9{r) obeying the boundary condition 9{R) = decays 
like 



= Cln- (2D) 
r 



(3D) 



(43a) 
(43b) 



as r — * R. 

The separation into two zones is also valid for solu- 
tions of the kinetic equation. It turns out that in the 
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run-out zone the space gradients are small and the diffu- 
sion asymptote (j3q) with 8 given by ( |43|) is reached. 

In the centre of the sample l n does not depend on n, 
because the solution is azimuthally symmetric. Taking 
into account that l\ is an odd function of n and using the 
constraint (40), we can define the value of l n at r = 
with a single parameter 8q- 

h(0)=0, / 2 (0) = isinh6>o, l 3 (0) = cosh0 o . (44) 

In the next subsection we find the exact solution of the 
kinetic equation for 8q = ln(l/7)- Although it does not 
match the asymptote (|38|), fl43| ) at large r, it plays an im- 
portant role. Any solution of the kinetic equation that 
starts at 8 < 111(1/7) eventually approaches the diffu- 
sion asymptote, while the one that starts at 8 Q > 111(1/7) 
does not. Therefore this exact solution is a separatrix. 

Now we are ready to describe the shape of the solution 
that obeys the boundary conditions (E3) and reaches the 
diffusion asymptote at large r. If 7 <C 1, it starts at 80, 
being only slightly smaller than lnl/7, and, therefore, 
runs close to the separatrix up to the large radius r* 3> I. 
For even larger r the deviation becomes significant and 
our solution crosses over to the diffusion asymptote j38| ) 
with 8 given by (ff3|) . The cross-over region is of the order 
of the mean free path I and is much smaller than both 
the reaction and run-out zones. We match the separatrix 
in the reaction zone with the diffusion asymptote in the 
run-out zone keeping the mean value (I) continuous at 
the point r* (see Fig. ||), thus finding both the position 
of the cross-over r* and the coefficient C in Eq. (43). 



B. Solutions 

An azimuthally symmetric solution of the kinetic equa- 
tion depends only on the radius r and the angle <p between 
the radius and the vector n. Therefore, the space deriva- 
tive in the kinetic equation has the form 

d sine/) d ( d , 

(45) 



nV 



cos< 



dr 



ds 



where the impact parameter p = r sin <fi and the distance 
along a straight line trajectory s = r cos <p are introduced 
(see Fig. |). 

If the parameterisation (|32| ) for the matrix g n is used 
together with the conic basis for the vector l n (see ap- 
pendix |a|) : 

h = k\i h — ih = k+ + "f, I3 + ih = k- + 7, (46) 

the kinetic equation is simplified to the form (in this sec- 
tion we measure all distances in the units of the mean 
free path I): 

'' ' =-(k+)h, (47a) 



Os 
d 



ds 

r„2 



k- 



(fc-)fcl, 

k\ + (k + +7)(fe_ +7) = 1. 



(47b) 
(47c) 



The condition at the origin has the form: 

fci (0) = 0, fc+ (0) = e e ° - 7, k_ (0) = e- 00 - 7. (48) 

The separatrix solution sta rts a t 80 — 111(1/7) an d there- 
fore, k-(0) = 0. Equation (47b) for fc_ now gives fc_ = 
for all r. Taking into account the strong inequality 7 <C 1, 
we simplify the system (EtI): 



dk+ 



ds 

k\ + 7 fc+ = 1 



(k+)k u 



(49a) 
(49b) 



and obtain a closed integro-diffcrcntial equation for k 
dk 



ds 



= -(fc + ) v /l- 7 /c+ 



The solution of this equation is given in appendix 
distances r > 1 it has the form: 



7*l, 



a 2 , , 6 ab (f — 101) In cot f 
— ln cot ■ " 



a = 2 < ln cot — 
2 



2 

-1/2 



sin ( 



4 f 1 (2D) 
n{V3 (3D) 



(f - |(/.|) ln cot § 



Sill ( 



4 In 2 



r 
95 



(2D) 
(3D) 



(50) 
|; at 

(51a) 
(51b) 
(51c) 



The divergence of k+ at <f> = is cutoff at the the angles 
(j> ~ 1/r 2 . Since all integrals with k + converge, the exact 
form of this cut-off is relevant only for the criterium of 
validity of the ballistic treatment (see section VII). The 
firstot term in ( |51a| ) has zero average and does not de- 
pend on r, while the mean value (fc + ) = a/(jr) decays 
when r — > 00. Using Eq. (49b), we find the expressions 
for the component k\ and the current J = (cos^fci) 



fci = ln I cot — I 

2 1 2 1 

a I 6\ 1(2 (2D) 

J =2\ COsWc0t 2 l / = n\V3 (3D) 



(52a) 
(52b) 



Note, that the the current J don't depend on the coor- 
dinates. 

We are looking at a solution of Eq. ( f47| ) that transfers 
to the diffusion regime at a certain radius r* 3> 1. The 
transition region has a width of the order of unity, and, 
therefore, both the mean value (&+) and the current J 
change negligiblely across the cross-over region. Com- 
bining Eqs. (H), (||) and © we obtain in the diffusion 
region: 



(*4 



a 

27 



(2D) 



C T CI .„„. 

= exp— , J = —n (3D). 



(53a) 
(53b) 



G 



Using continuity of J and at the point r* we find 
the values of the parameters C. r ■ • 



^1 In 1 

- 4 In f 



in ■ 



3^3 



7' 



c 



C = -4= In 2 i 

3\/3 7 



(2D), 
(3D) 



(54) 



The sketch of the solution is given on the Fig. 



C. Long-time Asymptote of Conductance 

We begin with the self-consistency condition ( |35| ) 
the conic coordinates (^) it looks as follows: 



In 



At 
irk 



dr (k + - k-) + 2-y 
V 2 



dr (k + 



(55) 



where it is taken into account that for solution (p 



the inequality fc_ <C k + holds. The last integral can 



be calculated using the continuity condition (41), where 



(k + — fc_)/2 can be replaced by k + /2. Integrating 
Eq. (41), we express the integral (5^) through the total 
current J through the outer boundary of the sample: 



J(r)dS 



r=R 



21 



dr(k^ 



(56) 



Although the integral (]55|) contains comparable contri- 
butions from both the reaction and run-out zones, with 
the help of Eq. ( |56| ) it can be expressed through the cur- 
rent in the run-out zone only. Using the asymptote (IBT 
we finally arrive at the expression for 7: 



7 = — pp/- In f 

' 2 yt t \pfIt 



7 



2T ln 2 



t 



T(p F lf 



(2D), 
(3D) 



(57a) 
(57b) 



The long time asymptote of the conductance G(t) i s dc - 
termined by the value of the action F from Eq. (151) 
for solution ( |5l| ) - ( |34|). This action is a sum of two 
contributions: from the run-out zone (-Fd) and from the 
reaction zone {Fb). The action is dominated by Fd' 



ttvD 



D 



dr 



ti 2 vDC 2 



ln(R/r* 
2/n 



(2D) 
(3D) 



(58) 



Substituting the values of the parameters ( |54| ) we obtain 
Tig In 2 tj (rg 



Fn = 



In R/l 



(2D) 



Fn = 



9^3 



(p F l)*\n 3 — (3D). 



(59) 
(60) 



where g — 2T:fivD is the dimensionless conductance. The 
contribution from t he re action zone Fb is calculated, us- 
ing ballistic action ( |l5b| ), in appendix ^: 



F b ~(j> F lf\n 2 



ln(i?/0 
t 



( PF l) 2 T 



(3D) 



(61a) 
(61b) 



Comparing Eqs. fl59|), ( |60[ ) and (^) we see that the ac- 
tion is dominated by the contribution from the diffusion 
zone. In two-dimensions it comes from the whole diffu- 
sion zone -c r -C R because the integral in Eq. j5^ ) 
diverges logarithmically. In three-dimensions only the re- 
gion near r = of width of the order of is important. 
Since 3> I, the contribution from the crossover region, 
which has the width of the order of I, can be neglected. 
The long time conductance asymptote is given by 

G(t)=exp{-F D ) (62) 

with the action Fn defined Eqs. (|59|) and m 



V. LONG-TIME ASYMPTOTE FOR 
CONDUCTANCE IN DISORDERED WIRE 



In this section we consider a ID thick wire of length L 
and cross-section w ( w <C I <S L) with specular bound- 
ary conditions and assume that the distribution function 
g n is uniform across the wire. For not very long times 
the diffusion equation (|37|) is valid and has the solution 



ifl 



£ = £/logiA. 



(63) 



The space gradient V9 = l/£ does not depend on the 
coordinate x and is smaller than l/l as long as the time 
is shorter than HA -1 exp(L/Z). For longer times the dif- 
fusion regime breaks down simultaneously in the whole 
wire. The separation into reaction and run-out zone is 
still valid, and the term in u> in the kinetic equation is 
important in the reaction zone near the centre of the wire 
and can be neglected elsewhere. Unlike the 2D and 3D 
cases however, the space gradients in the run-out zone do 
not decay towards the outer ends of the wire. We now 
present a solution of the kinetic equation in the run-out 
zone which is valid for arbitrary gradients. 

The distribution function g n depends on the coordi- 
nate x along the wire and the angle <p between the axis 
of the wire and the direction n of the momentum (see 
Fig. |l|). In these coordinates the kinetic equation (||) 
takes the form: 



2lco8^ = [(7A-(ff)) )fl ] 



(64) 



In the run-out zone the term in 7 is negligible and, using 
the conic coordinates, from appendix ^ as 7 — * 0, we 
rewrite Eq. (p3): 



7 



lcoscj)—k ± = T{k±)k x , 



k + k- 



1. 



(65a) 
(65b) 



Dividing both sides of Eq. (|65a| ) by cos <p and averaging 
over cj), we find a closed equation for (fc±) with the solu- 
tions: 

(fc±) =qexp( T 0), (66) 
where q is a constant and 9 obeys the equation: 
d9 . fci . 



dx 



cos< 



(67) 



In one-dimension the current conservation law ( |42| ) re- 
duces to J = (ki cos 4>) — const and suggests that k\ 
does not depend on x. Therefore, Eq. (K37n gives 



I— — = const = 9 
dx 



(68) 



and, substituting this back into Eq. (65a), we obtain the 
solution in the factorised form: 

k± = P(4>) (k± ) , fci = 6'f3{<j)) cos 0, (69) 

where the function (3{4>) is determined from Eq. ( |65h| ) 

P=[q 2 + {6') 2 cos 2 (70) 

Finally, the constants 6*' and g are related by the condi- 
tion 



1 f 2 " 

1 = w = ^i 



vV + {9') 2 cos 2 



(71) 



Using the asymptotic values of the elliptic integral in 
Eq. ([71]), we find 



<1 = 



40'exp( - 2|1 



9'< 1, 
9'> 1 



(72) 



and obtain expressions for the current J = (ki cos < 



J = 



f {l-47r|0'|exp(-7r|0'|)} 



9' < 1, 
9' > 1 



The above solution is not valid in the reaction zone 
in the vicinity of x = 0, whose contribution to both the 
action and self-consistency condition is negligible. The 
boundary condition for the ballistic problem is deter- 
mined by the requirement that the distribution function 
matches the one in the bulk electrodes where g n = A. 
For small gradients (6' -^i 1) the angular dependence 
of the distribution function in the wire is weak and the 
above condition is fulfilled provided 9(±L/2) = 0. On 
the other hand, in the ballistic regime (Q 1 3> 1) the dis- 
tribution function in the wire strongly depends on the 



angle <j) (see Eq. ([70|)). Thus, there is a cross-over region 
near the outer ends of the wire where the solution of the 
kinetic equation deviates from Eqs. (|66|), ( |69| ) and (|7C|). 
We consider the wires which are long enough to neglect 
the change of 9 in the cross-over region. Therefore, the 
boundary condition 9(±L/2) = is still valid in the bal- 
listic regime, and the function 9{x) can be presented in 
the form (||) with £ = 1/9' and 9 = 9'L/(2l). 



A. Action and self-consistency condition 



The self-consistency condition ( fig ) gives 



t A L 1 f\9'\L 
57,exp 



nh I \e> 



21 



(74) 



and should be compared with the continuity equation 
( [56| ) which in the ID case reads: 



J 



2ik exp 



\0'\L 
21 



(75) 



Using Eq. @ for the current and Eqs. Q, and @ we 
express 7 and 9' through t: 



\9'\ 



21 



L 

2-kTi I 



1 = 1Kl j= t 



(l n ^ + l n l n !^y ( 76a ) 
rf2<7lni ^«exp(f), 



2N 



H» ex p(f ' 



(76b) 



where g = 2-nhvDw I L is dimensionless conductance and 
N = wpp I nh is the number of transverse channels, which 
is equal to ballistic conductance. The calculation of the 
W-term in action (15b) for the solution ( |69| ) is performed 
in appendix ^ and gives: 



9' J 



(73) Evaluating the other terms in the action we get: 



F 



2A- 



-(2J9' + q 2 -l) 



(77) 



(78) 



Thus, the long time asymptote of the conductance is 
given by: 



eL/! « « NeL/l 



G(t)^exp{-.gln 2 (^.)} 1 « ^ « e&' 
G(t) - (H) 

The last inequ 

semiclassical approximation (see section VII) 



(79) 



The last inequality in Eq. d79) ensures the validity of 

an (see 
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VI. PHYSICAL PICTURE OF TRAPPING 

In this section we consider the time dispersion of the 
conductance G(t) in a purely one-dimensional wire. This 
problem was solved by Altshuler and Prigodin jl2| who 
obtained the long time asymptote in the form 



G{t) ~ exp ( --J-ln 2 tA 



(80) 



which can be treated as a limiting case of a multi-channel 
formula for a thick wire (see table |lj), assuming that in 
a one-channel case g is given by l/L. Formula (|o|) can 
be understood as a probability of an optimal potential 
fluctuation that traps an electron of Fermi energy Et for 
time t. In a weak potential U(x) <C Ef the wave function 
can be presented in the form 



#(3) = (t>+(x)e tpFX + 4>^{x)e 



(81) 



with the amplitudes <f>±(x) changing slowly: V</>± -C 
Pf4>- Let us consider a quasi-stationary state obeying 
the open boundary conditions 



4>+(0) = (j)-{L) = 0, 



(82) 



which correspond to the outward flow of current through 
the ends of the wire. The life time of such a state is 
inversely proportional to the outward current 



t 



fdx\*\ 



^(|^_(O)| 2 + |0 + (L)P)- 



(83) 



The maximum delay time is achieved when the currents 
through both ends are equal (</>-(— L/2) — <f> + (L/2)). 
Fixing normalisation by 



we reduce Eq. 



IMo)| = \ML)\ = i 



to the from 



(84) 



(85) 



which resembles the self-consistency condition (Q) ob- 
tained for arbitrary dimensions. 

To achieve life times t > li/A the wave function must 
grow towards the middle of the wire; assuming the growth 
is exponential, \& ~ exp[(L/2 — |x|)/£], we obtain for the 
localisation length of the quasi-stationary state 



L 



lniA 



(86) 



A typical random potential U causes one-dimensional 
wave functions to be localised with £ ~ I. The shorter 
localisation length £ <C I corresponds to life times longer 
than ?iA _1 exp(L/7) and can be achieved in the potential 



with the additional 2pp-Fourier component having the 
amplitude 



2hv f 



The probability of this potential realization in given by 
the Gaussian distribution: 



exp - 



■KVT 



J U{x) 2 dx^j ~ exp (-^ In 2 tA 



(89) 



and coincides with ( |80|) with a correct numerical factor in 
the exponent. We therefore conclude that the states with 
long life times are locked by the Bragg reflection and can 
be found with the probability proportional to that of po- 
tential fluctuation with the Bragg mirror of appropriate 
strength. 

We believe that the same mechanism is responsible 
for nearly localised states in multi-channel wires and the 
samples of higher dimensionals. In multi-channal cases, 
however, adding a single 2p^-Fourier harmonics cannot 
localise the wave function, because the random part U 
mixes different directions of the momentum. To localise 
the wave function in a two- or three- dimensional sam- 
ple, the potential fluctuation should be effective for all 
directions of the momentum; so we expect it to have the 
form 

U(x) = J dfl n U n (r) cos(2p F nr) (90) 
with the amplitude U n (r) slowly depending on n and r. 

VII. DISCUSSION 

This paper continues our study of the conductance 
G(t) asymptote at long times t that began in Ref. [0. 
We use the steepest descent approach which enables us 
to obtain g(t) for different ranges of time. The purpose 
of this section is dual: we want to analyse the restrictions 
of our treatment and compare our results with those in 
the literature. 

Let us first discuss the conductance G(t) of a thick wire 
made from a two-dimensional strip of length L 3> I and 
width w < I <C L. The solution g(x, <f>) of the kinetic 
equation has typical gradients (see Eq. (|77|)): 



dg & \n{tA/K) 
dx I L 



(91) 



At long times the distribution function changes rapidly 
with the angle <f> acquiring a sharp maximum at <p = tt/2 
with the width 5<f> - exp(-7r|0'|/4) (see Eq. @-@). 
The semi-classical approximation remains valid if this 
width is much larger than the diffraction angle 5<j)Q ~ 
Pfw. This imposes an upper limit on the gradients: 



U(x) = U + U cos(2p F x). 



(87) 



6' < \n(p F w/h) 



(92) 
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Thus, the ballistic asymptote is valid only for times 
smaller than 



t 



Q 



A~ 1 p F we L/l 



(93) 



Since the width of the wire is limited by the condition Jig] 
w < I the gradient 9' at times t ~ tQ is still much smaller 
than pp. Thus, an intermediate interval of gradients 
arises 

\n(lp F ) < 6' < l PF (94) 
which corresponds to the interval of delay times: 

~]n{«ln(fA)«~. (95) 
/A A 

We assume that the time dispersion of the conductance 
in this region can be recovered by some kind of a ballistic 
treatment with diffrction proper accounted for. 

A similar phenomenon restricts the range of applicabil- 
ity of ballistic asymptotes in two- and three- dimensional 
samples. Analogously to the one-dimensional case, the 
distribution function g(4>,r) in the reaction zone has a 
sharp feature at <fi — with the width 8(f) = l/r. Diffrac- 
tion can be neglected when 8(f) ^> X/l. Therefore, the 
reaction zone radius r* must obey the constraint: 



r, < l 2 /\. 



(96) 



Substituting the radius of reaction zone from Eq. (^4|) , we 
obtain the upper time limit t for which the semi-classical 
approach is still valid: 



t«.t Q ~ hA 1 exp 



n 



l/R 
(l/Rf 



(2D) 
(3D) 



(97) 



Similarly to the one-dimensional case the gradients of g 
at the time tQ are still much smaller than p F , being only 
of the order of l/l. The liming times tQ for a different 
geometries are summarised in table |. 

As it has been mentioned in the introduction, the 
whole field was pioneered by AKL. They added high pow- 
ers of frequency and gradients of Q to the diffusive a- 
model and analysed the renormalization flow of the cor- 
responding coupling constants in two dimensions. This 
gave the growth rate of the coefficients C n in the expan- 
sion 



G(w) 



C n u n . 



(98) 



The Fourier transform of Eq. ( |98|) with the AKL asymp- 
tote for C n leads to the result presented in Table |[ AKL 
also put forward a general conjecture that the logarith- 
mically normal asymptote was valid for all dimensions. 

As one can see from table [| there is a variety of dif- 
ferent regimes in the time dispersion of the conductance. 
The AKL procedure predicts one of them and fails to de- 
scribe the others. It is also difficult, using this procedure, 
to find out the criteria when this or that results are valid. 



Another attack on the problem has been recently car- 
ried out by Mirlin ||. He used our optimal fluctuation 
method for the diffusive er-model in dimensions two and 
three and supplied Eq. (|^) by somewhat arbitrary the 
conditions near the origin. In this way he obtained the 
asymptotes of conductance for 2D and 3D samples. The 
reason for his success lies in the fact established in this 
paper that the action for the dimensions d > 2 is dom- 
inated by the contribution from the run-out zone. The 
solution of the ballistic problem confirms the result of 
Ref. m for two-dimensional conductance. In the three- 
dimensional case we obtain the numerical coefficient A in 
the expression G(t) ~ exp(— A\v? t). Although the form 
of the density-matrix fluctuation suggested by Mirlin de- 
viates strongly in the reaction zone from the correct one, 
the method of his work could be useful for qualitative 
estimates. 

It should be noted, however, that the ballistic treat- 
ment is needed to find the upper limit tQ on the delay 
time interval where these results are valid. We regard the 
occurance of the the new range of times t > tQ where the 
diffraction effects are important we regard as one of the 
most interesting results of this paper. 
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APPENDIX A: CONIC COORDINATES 

Throughout this paper we use a convenient parametri- 
sation of the bosonic sector of the distribution function 
g n . To obtain it, we introduce a basis in the space of 
traceless 2x2 real matrices: 



ei = oi, e± 



o- 3 ± ia 2 



(Al) 



The new unit vectors have the following properties: 



e\ = 1, e± = 0, [e±,ei] = ±2e±, 
[e+,e_] = -ei, {e + ,e_} = 1. 



(A2a) 
(A2b) 

Now, instead of parametrisation ( |32"| ) for the g^, we use: 

9n = ® (Mn)ei + h{n)e+ + /_(n)e_) + 

H — - ® (Zi(-n)ei + l+(-n)i + + Z_(-n)e_) , (A3) 



where, due to requirement (p9|), all components l\ and l± 
are real. The kinetic equation can be further simplified 
by introducing the functions 
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k% = h; k + = l + -j; k- = I- — 7 (A4) 
The constraint g 2 — 1 now can be written as: 

fc? + (fc++7)(fc-+7) = l, (A5) 

and the kinetic equation ^ in coordinates ( [45| ) has the 
form: 



d_ 

da 

d_ 

ds 
dki 



k + = -{k + )k u (A6a) 
k- = (fc-)fci, (A6b) 
= -k+(k-) + <fc+)fc- + j(k+ - k-). (A6c) 



It is more convenient to use the constraint (A5) instead 



of the last equation in this set (see Eq. (47) in the text 



APPENDIX C: CALCULATION OF BALLISTIC 
ACTION IN THE REACTION ZONE 



The contribution from the reaction zone to the action 



is given by Eq. (15t) and contains two terms: 



jr 1= ™* Wl F 2 



8r 



<5> 2 ; 



(Cl) 



the term with to dissapears after integration over oj in 
Eq. ©. 

To caluclate the W term we use Eq. ([l4]) and present 
the solution of kinetic equation in the form g n = U AC/ -1 . 
Using parametrisation ( |32| ) we find for the bosonic block 
of g-matrix: 



APPENDIX B: INTEGRO-DIFFERENTIAL 
EQUATION FOR SEPARATRIX 

Solving Eq. j50| ) with respect to fc+ and considering 
its average value (fc-|_) as a given function of radius r, we 
obtain: 



1 - 7 fc+ = 



7 



(k+)ds' 



(Bl) 



By taking the average of both the right- and left-hand 
sides, this expression is reduced to an integral equation 
for A(r) = 7(fc+): 



1-A(r) = 7 



where 



r cos <fi 



ds\ 



r 2 sin 2 1 



(B2) 



A T27T d4 

\><t> JO 271- 
/\ r 71 " sin d>dt 

00 = Jo ~2 



(2D) 
(3D) 



The solution of Eqs. ( |B2[ ) has the following asymptotes 
A(0) = 1; 



A(r) = ^ + A r »i j (B3 ) 



where the constants a and b are given by Eqs. (|5l|) in the 
text. 

If r sin <p ^> 1, then the asymptote of ( [B3| ) can be sub- 
stituted into kernel of Eq. (Bl), and give the result 



7*l, 



° 2 1 2 , 
1 In cot • 



a6 ln (| -|^|)| cot I 



2r 



sm i 



(B4) 



Expression (B4) is valid only if r\ sin0| 3> 1 and is free 
of singularities in this region. 



9 = 



U^azU- 1 ^) 

-U{-K-4>)a z U- 1 {-K 



where 2x2 matrix U is given by: 



/ = cxp j --^e+^J ex P (~2~ <Tl 1 ' 



(C2) 



(C3) 



Using 



9 TJ ldki - TT 

7T U = TT e + U 

or 7 or 



we obtain from Eq. (JT3) : 



IV = s I dr(n^-) = 8 / <1SJ. 



(C4) 



where the integral in the last expression is taken over the 
boundary of the reaction zone. This gives 



ln(£/£T 
y\n(RJl) 



(p F l) 2 \n 21 * 



(2D) 



(C5) 



which is one power of the logarithm smaller than the 
contribution from the run-out zone. 

The calculation of the term F2 is strait-forward: 



F 2 = 



8r 



dr- 



(^) 2 m((^) (3D) 

(C6) 
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APPENDIX D: CALCULATION OF W TERM 
FOR ONE-DIMENSIONAL SOLUTION. 

To employ the expression ( |l4| ) for the W-term we 
should find a decomposition g n = UAU^ 1 . Combining 
Eqs. ( |66| ) and (|6|) with the parametrisation (|?2|) and for- 
mulae from app end ix |a|, we write the bosonic block of g n 
in the form of (|C2D, where the 2x2 matrix U is given 
by: 



U = cxp ( ^c^ J exp (-*|cr y ) , 



where x — arcsin (Q' f3(<j)) cos (f>) does not depend on x. 
The straight-forward application of Eq. ( |l4|) gives: 



W = -8W 



9V 
I 



(D2) 



where the minus sign comes from the super-trace of the 
bosonic block. The result of averaging over n is expressed 
through the current J = (cos^fci). 
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t D < t < h/A 


h/A < t < t b 


t b A 


t b < t < f Q 


t Q A 


lD a 


t/t D 


g\n\tA) 


exp(L/Z) 


2 AH n i A 


Nexp(L/l) 


2D b 


irt/t D 


Ag ln(tA) 


{L/lf 


TJ ln 2 (t/ S r) 
71 9 2\n{R/l) 


ey:p{p F l/h)ji 


3D C 


llt/t D 


none 


1 




t 




exp(p F l/K) (j^ 3 







TABLE I. Long time asymptotes of the conductance G(t). The cells contain formulae for — lnG at different time intervals 
for one- two- and three- dimesional samples. The diffusion time to is defined as L 2 /D (Id) and R 2 /D (2d, 3d); the diffusion 
coefficient D is given by \v\t (Id, 2d) and \v"pT (3d); r is the mean free time and A = (vV) -1 is the mean level spacing. 
a strip of width w and lenth L, N = wpF/(^h) is the number of ballisitc channels 
b disk of radius R 
c ball of radius R 
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FIG. 1. The geometry of the Id sample (a) and 2d and 3d 
samples (b) 

FIG. 2. Sketch of kinetic equation solutions in two- and 
three- dimensional samples, (dotted line is the separatrix so- 
lution) 
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FIG. 1. The geometry of the Id sample (a) and 2d and 3d samples (b) 




FIG. 2. Sketch of kinetic equation solutions in two- and three- dimensional samples, (dotted line is the separatrix solution) 



